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ABSTRACT 



The numerous and massive young star clusters in blue compact galaxies (BCGs) 
are used to investigate the properties of their hosts. We test whether BCGs follow 
claimed relations between cluster populations and their hosts, such as the the fraction 
of the total luminosity contributed by the clusters as function of the mean star forma- 
tion rate density; the V band luminosity of the brightest youngest cluster as related to 
the mean host star formation rate; and the cluster formation efficiency (i.e., the frac- 
tion of star formation happening in star clusters) versus the density of the SFR. We 
find that BCGs follow the trends, supporting a scenario where cluster formation and 
environmental properties of the host are correlated. They occupy, in all the diagrams, 
the regions of higher SFRs, as expected by the extreme nature of the starbursts oper- 
ating in these systems. We find that the star clusters contribute almost to the 20 % of 
the UV luminosity of the hosts. We suggest that the BCG starburst environment has 
most likely favoured the compression and collapse of the giant molecular clouds, en- 
hancing the local star formation efficiency, so that massive clusters have been formed. 
The estimated cluster formation efficiency supports this scenario. BCGs have a cluster 
formation efficiency comparable to luminous IR galaxies and spiral starburst nuclei 
(the averaged value is ~ 35 %) which is much higher than the 8-10 % reported for 
quiescent spirals and dwarf star-forming galaxies. 

Key words: galaxies: starburst - galaxies: star clusters: general - galaxies: irregular 
- galaxies:star formation 



1 INTRODUCTION 

Luminous blue compact galaxies (BCGs) are star-forming 
systems with high specific st ar formation rate (|Ostlin et al.l 
I2OOII : ISoUerman et al.l 120051 ). It has been suggested that 
BCGs may have accounted for ~ 40% of the total star for- 
mation rate (SFR) density at redshift 0.4 <z< 1.0 play- 
ing an importan t role in the star for mation history (SFH) 
of the Universe (|Guzman et al.lll997h . However, in the lo- 
cal Universe, their contribution has dropped drastically and 
luminous BCGs with high SFR have b ecome rare objects 
l|Guzman et all Il997l : IWerk et allbOOll ). At a distance of 
100 Mpc (z-^O.OS), we count only a handful number of BCGs 
with a SFR higher or close to 5 M^yr"^. Among the s ystems 
included i n this work, Haro 11 (JAdamo et aklbOlOaTl. ESQ 
185-IG 13 (|Adamo et aLlbOllal ). and Mrk930 (JAdamo et al.l 
uOllbl ) are starburst BCGs with SFRs exceptionally high 
for local irregular galaxies. Local luminous BCGs display 
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physical conditions (low metallicity and dust content), mor- 
phologies (compactness of the starburst re gions), and feed- 
back mechanisms (e.g. lOstlin et al.l [2009|) similar to their 
high redshift counterparts. 

Numerical simulations and theoretical arguments based 
on the Lambda Cold Dark Matter model of the Universe 
predict that smaller galaxies formed first and were then 
accreted into more massive systems (hierarchical growth , 
ICole et al.|[2OO0l : ISpringel et aUbOOSl : iPapovich etlLlbOOSl ) . 
Therefore, it is expected that the primordial galaxies were 
chemically unevolved dwarf galaxies, the so-called "building- 
blocks" systems, which accrete and merge into more massive 
units. Recent observations of galaxies at very high redshifts 
(z~ 7, i.e. at the reionization epoch) have revealed compact 
single or double nuclei s ystems with extended nebular fea- 
tures (jOesch et al.ll2010l ). Some of these z~ 7 objects have 
also been detected at longer IR wavelengths. These data 
have been used to constrain and study the spectral energy 
distributions (SEDs) of these high -z systems. The infe rred 
SFRs are between 5 to 20 Moyr"^ (JGonzalez et al.ll2010l ). In 
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general, they have estimated stellar masses of 10*~^ Mq and 
fainter UV luminosities than l ower redshift Lyman break 
galaxies (JFinkelstein et al.ll201Ct l. i.e. their properties resem- 
ble those of the luminous BCGs. 

The studies of BCGs can, therefore, give important 
insights in understanding how star formation proceeds in 
dwarf starburst galaxies and possibly in high-redshift sys- 
tems, with spatial and spectral resolution impossible to 
achieve for the latter with the current facilities. High- 
resolution imaging data of BCGs revealed that the starburst 
regions in these galaxies ar e formed b y massive a nd young 
star cluster co mplexes (e.g. lOstlin et a l. 2003,; Ad amo et al.l 
l2010al . l2011al lbh. The peak age of the star cluster distribu- 
tions and the estimated ages of the starbursts in these sys- 
tems are in good agreement. The host morphologies suggest 
that the galaxies have recently undergone a merger or inter- 
action event, which has likely refurbished the galaxy with 
metal-poor gas, and triggered a vigorous starburst episode. 

Star clu sters are a natural ou tcome of the star formation 
process (see lLada fc Ladall2003l for a review). However, for 
very young stellar systems (e.g. a few Myr or less), it is not 
trivia l to make a clear definition of a cluster (JBressert et al.l 
|2010| '1. Usually, it is assumed that stars form in a clustered 
fashion and that, after roughly 1 Myr, 90 % of clusters are 
destroyed due to gas expulsion. iBressert et al.l (|2010|) sug- 
gest that all the stars form in a continuos hierarchy, even 
the dense regions and estimate that, in the solar neighbour- 
hood, only ~ 26% of young stellar objects (new born stars) 
are located in denser regions, i.e. embedded star clusters. 
The remaining stars form in associations and agglomerates 
following a hierarchical contin uum distribution with c l uster s 
at t he bottom of the process. IPortegies Zwart et al.l (|201Cl ) 



and iGieles fc Portegies ZwartI (2010|) defined an empirical 
relation to separate clusters from loose associations. For the 
latter the crossing time is much larger than the age of the 
stars, i.e. they are dynamically unbound systems. They ob- 
served that the division between bound and unbound sys- 
tems becomes more clear once the gas-expulsion phase is 
over (~ 10-20 Myr). 

Cluster formation at high-redshift is almost unknown. 
In the local Universe, we observe the evolved counterparts, 
i.e., the globular clusters (GCs). However, there is no di- 
rect evidence that the young star clusters we observe lo- 
cally will eventually evolve into GCs (|de Griis fc Parmentierl 
[20071). If one looks at the number of GCs per mass (lu- 
minosity) bins (i.e., dN{M)/dM = Ch'^) the distribution 
can be fitted by two power-laws with a slope of ~ —2.0 
at the massive (luminous) ends, a break corresponding to 
the M c ~ 2 X 10^ M0, and a flattening (77 ~ -0.2 which 
variates from system to system) at lower masses. Sev- 
eral possible scenarios are addressed by theoretical studies 
to explain the origin of the GC mass (luminosity) func- 
tion: dynamical evolution, i.e. due to preferential disrup- 
tion of the low mass systems; or a primordial origin, i.e. 
the GC mass function has been established at th e time 
when the GCs formed (|de Griis fc Parmentierl 120071 ). Cos- 
mological simulations seem to su ggest that GCs may have 
formed in dark matter ha los (e.g. iKravtsov fc Gnedinll2005l : 
iMashchenko fc Silld |2005| V If this is the case, there is no 
connection between the YSCs forming at the present time 
and the ancient GCs. However, some recent works (see 
iBrodie fc Straden 120061 . for a review on GCs) suggest that 



cluster formation in dwarf galaxies may be the key to ex- 
plain observed properties of the ancie nt GCs (blue versus 
red, i. e. metal poor versus metal rich) . iMuratov fc GnedinI 
(|2010h recovered a bimodal metallicity distribution as a 
product of cluster formation in different phases of galaxy 
evolution. Cluster formation in dwarf galaxies produced the 
blue, metal-poor GCs. Dwarf systems were successively ac- 
creted to form more massive systems. During the merger and 
formation of these massive galaxies the more metal-rich clus- 
ters were formed. Accretion of dwarf galaxies together with 
their GC systems i s also one of the proposed scenarios by 
IChies-Santos et al.l (|201l|) to explain why the younger GCs 
in SO type of galaxies appear blue instead of t he expected 
metal-rich populations (jBrodie fc Stradeiil2006l ). 

In the present work, we will investigate how cluster for- 
mation has proceeded in BCGs. We will test whether BCGs 
follow the cluster-host relations available in the literature 
and constrained using local star-forming galaxies, like spirals 
and dwarfs. These results will be used to constrain the envi- 
ronmental properties of BCGs and whether star formation 
operates on similar modes even under extreme conditions. 

The paper is organized as follow: In Section 2, we 
present the 5 BCG targets used in this work. In Section 
3, we first discuss the uncertainties which affect the anal- 
ysis. In the second part of this section, we show the three 
cluster-host relations including the BCG sample. A discus- 
sion of the results is presented in Section 4. Here, we also 
discuss possible similarities between BCGs and high redshift 
galaxies. Conclusions are summarized in the last section. 



2 THE DATA 

In this section we shortly introduce the BCGs included in the 
analysis. Some of these BCGs have been stud ied in a series 
of 3 pa pers: Haro 11 analysis is pres ented in Adamo et al.l 
(|2010al ): ESQ 1 85-IG13 (ESQ 185) inlAdamo et all (|2011al '): 
and Mrk930 in lAdamo et all (|2011bl '). We refer to those pa- 
pers for details on the analysis of the data used in this work. 



2.1 ESO 338-IG04 

The analysis of the star cluster pop ulation in ESO 338 - 
IG04 (ESO 338) has been presented in lOstlin etHI (|2003h . 
The masses have been obtained from models that as sume 
a Salpeter initial mass function flMF. ISalpeteilll955l ). We 
show, in Figure [T] the cluster formation history during the 
last 40 Myr of galaxy evolu t ion. U sing the age and mass 
estimates from lOstlin et al.l (|2003l ). we assumed that the 
analysis is complete in detecting clusters more massive than 
5 X 10"^ Mq formed during the last 40 Myr. A power law 
cluster mass function with index —2.0 has then been used 
to extrapolate the total fraction of mass in clusters includ- 
ing obj ects with mass e s betw een 10^ <M< 5 x 10^ M©. Fol- 
lowing |Weidner_et_alJ (|200J), we assume for simplicity that 
a cluster population forms every 10 Myr and estimate the 
clust er formation rate (CFR) in the galaxy. In agreement 
wit h i Ostlin et al.l (|2003h . we observe a cluster formation en- 
hancement between 20 and 30 Myr ago. However, the peak 
of cluster formation is younger than 10 Myr, similarly to 
what is found in Haro 11, ESO 185, and Mrk930. ESO 338 
is forming roughly 1.6 A/0yr~^ of stars in clusters. 
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Figure 1. Cluster formation rates during the last 40 Myr of star- 
burst activity in ESO 338. The filled dots connected by the green 
dotted line show the observed CFR derived from clusters more 
massive than 10^ Mq. The squares show the derived CFR, if the 
total mass in clusters less massive than 10'' Mq is extrapolated 
using a CMP with index — 2.0 down to IO^Mq. The data used 
are from lOstlin et al.l | |2003|) . 



We defined the cluster formation efficiency (CFE, indi- 
cated hereafter also as F) as that the fraction of stars formed 
in star clusters. If we compare the CFR to the mean SFR 
happeni ng in the galaxy , we obtain the CFE=CFR/SFR (see 
GIO and lBastianll2008l here after BOS). Using a SFR of 3.2 
Moyr-^ (Ostlinet al.ll200lD . we find that F ~50 ± 10 %, 
i.e. 50 % of the ESO 338 stars are formed in clusters. 
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2.2 SBS 0335-052E 

The nuclear region of this extreme metal-poor galaxy is 
dominated by 6 massive young star clusters which have 
been referred to in the literature as super star clus- 
ters (SSCs). T hese SSCs have been extensively stud- 
ied before (eg.. iReines et all l2008l : iThompson et all |2009| : 
lAdamo et al.ll2010bl among the most recent published works 
on this subject). 

In the present work we explore the lower mass clus- 
ter population of the galaxy. Part of this u nderlying cluster 
population has already been revealed by IPapaderos et al.l 
(| 19981 ). We have access to FVV (S'BC/F140LP from 
GO 9470, PI; Kunth) and optical images {ACS F220W, 
F330W, F435W, F550M from GO 10575, PI: Ostlin) for the 
galaxy. The reduc tion of the science frames is described in 
lOstlin et al.l (|2009l ) . The extraction of the cluster candidates 
has been done using the PyRAF package DAOFIND on the B 
band frame (F435W). The catalogue has been cleaned by 
eye of all the detections which didn't show a clear visual 
counterpart. With this first catalogue, we have done pho- 
tometry on all the frames from the FU V to the V (F550M) 
bands. Applying the same method (see lAdamo et al.ll2010al . 
l2011al lbl) we have included in the final catalogue only cluster 
candidates with detection in at least 3 filters and a photo- 
metric error (j„i < 0.2 mag. The photometric properties of 
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Figure 2. In the top panel the color-color diagram of the whole 
cluster population detected in SBS 0335-052E is shown. Blue dots 
indicate the six super massive star clusters, which dominate the 
starburst in the galaxy and have already been studied in several 
works (e.g. lAdamo et al.ll2010bl for references). The black dots 
represent the newly studied lower mass cluster po pulation. The 
evol utionar y track has a metallicity of Z = 0.0004 lAdamo et al.l 
("see l2010bl . for details). Some important ages are labelled. The 
extinction vector indicates in which direction the clusters move if 
a de-reddening of Ky = 1 mag is applied to the observed colors. 
A mean photometric error-bar is shown. In the bottom panel, we 
show the mass-age diagram of the star cluster population. The 
several lines indicate the detection limits reached in our observed 
data (shown in the inset). See main text for details. 



the final cluster population are shown in the color-color di- 
agram in Figure [21 The filter combination corresponds to a 
U-B (F330W-F435W) versus &B-V (F435W-F550M) 
color. The colors are not dereddened but clearly suggest that 
the cluster candidates are not older than ~ 100 Myr. 

The spectral en ergy distribut i on (SE P) fitting proce- 
dure is described in lAdamo et al.l (l2010al ). The used mod- 
els are presented in lAdamo et al.l (|2010bl ) . The output age 
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and mass of the clusters are shown in the mass-age diagram 
in Figure (2] In blue dots, we show the 6 SSCs previously 
analysed. The underlying black dots represent the low mass 
cluster population that is present in the galaxy, detectable 
at the detection limits imposed by the data. The masses 
are smaller than ~ 5 x 10* Mq and the ages are younger 
than 50 Myr. It is still under debate whether an old stellar 
population (older than 100 Myr) exists in SBS 0335-052E 
Ostlin fc Kunthll2001[l or if this galaxy has recently formed 



I 



Paoaderos et ah l lT^. In previous star cluster analyses of 



other BCGs (ESO 338, Haro 11, ESO 185, Mrk930), we have 
found a trace of some old GCs, supporting the evidence of 
an old underlying stellar population in those galaxies. In the 
case of SBS 0335-052E, however, the non detection of mas- 
sive GCs cannot prove/disprove either of the two proposed 
scenarios. We observe that, because of the detection limits 
(see inset in Figure [2]), our analysis is limited to GCs with 
masses higher than 5 x 10* Mq between 100 Myr and 1 Gyr 
and even more massive at older ages. Therefore, we cannot 
exclude that this galaxy has low-mass GCs. 

Using the same method as in the case of ESO 338, we 
infer the mass in clusters formed during the last 10 Myr 
in SBS 0335. We observe that the total mass contained in 
detected clusters younger than 10 Myr is 1.95 x 10^ Mq . The 
mass contained in the 5 SSCs (one of them is much older, 
e.g. from Ha equivalent width, the age is cons trained to ~ 13 
Myr is constrained, see lAdamo et al.ll2010bl ) is roughly 73 
% of this total mass (1.42 x 10^ Mq). Assuming a power law 
cluster mass function and that we are complete in detecting 
clusters more massive than 5 x 1O^M0, we estimate a total 
mass in clusters younger than 10 Myr of 6.4 x W^Mq. The 
observed CFR in systems more massive than M> 5 x 10^ Mq 
is 0.2 M0yr~\ while the extrapolated CFR (M> 10^ Mq) = 
0.64 Moyr"\ The F value of SBS 0335 is 49 ±12%, using 
a SFR of 1.3 MQ-yr"^ . These values have been estimated 
after that a correction for Salpeter IMF has been applied 
(see Section [3. l.ip . 



3 RELATIONS BETWEEN CLUSTER AND 
STAR FORMATION RATE 

The formation of a cluster appears to be correlated with the 
properties of the host galaxy. A common observation is that 
galaxy mergers produce more numerous and more massive 
clusters than quiescent spirals (|Larsenll2009l . and references 
therein). Sampling statistics, known also as size-of-sample ef- 
fect, (i.e., galaxies with a more numerous cluster population 
have higher chances to sample the cluster mass function, 
CMF, at higher masses) is a possib le explanation for this 
trend (|Whitmorell200d : lLarsenll2002l . hereafter L02). On the 
other hand, the host environment can pl ay its role in deter - 
mining the mass of the forming clusters (jGieles et al.ll2006r ). 
Numerical simulations of different host environments sug- 
gest that the shear in rotationally supported galaxies (i.e., 
spirals) acts on the collapse of the giant molecular clouds 
(GMCs), causing fragmentation and favouring the formation 
of the less clustered OB associations and low mass clusters 
IWeidner et al.ll20ld ). The lack of rotation in dwarf galaxies 
and high external pressures in merging systems favour the 
collapse of massive and gravitationall y bound cluster. Th ese 
two scenarios were also addressed bv lBillett et al.l (|2003) to 



understand cluster formation in dwarf galaxies. They stud- 
ied the star cluster populations of nearby dwarf galaxies 
and observed that not all the systems had bound and lu- 
minous clusters. However, some of them hosted one or a 
few very massive ones. They suggested that, with respect to 
spiral galaxies which form more clusters and can, therefore, 
sample the cluster mass function homogeneously up to high 
mass bins, the cluster formation in dwarf galaxies is pos- 
sibly dominated by stochasticity together with favourable 
physical conditions to form single massive clusters. 

Observed empirical relations between the properties of 
the young star clusters and the SFR in the host s u pport 
the size-of-sample effect scenario. iLarsen fc Richtlen (|2000l . 
hereafter LROO) first noticed that the fraction of luminosity 
contained in the young star clusters and the SFR of the host 
galaxy are correlated (Tl(U)-Esfr realtion), i.e., higher 
SFRs correspond to a more numerous cluster population 
(higher cluster formation efficiency). In a follow-up work, 
L02 found evidence of a positive correlation between the vi- 
sual luminosity of the brightest star cluster and the SFR in 
the host (MJ;"'^''*''''*-SFR relation). The relation between the 
two quantities can be understood if higher SFRs enabled the 
formation of more massive clusters. B08 enlarged the sam- 
ple of L02, including resolved close-by star-forming regions 
and luminous and ultra-luminous IR galaxies (LIRGs and 
ULIRGs). He noticed that the M^"«^''"=''-SFR relation holds 
over several orders of magnitude in SFR values, suggesting 
that the youngest brightest cluster is a fairly good indicator 
of the present SFR in the galaxy. These observed relations 
clearly point towards a scenario where the cluster formation 
is intimately correlated with the star formation process, or 
in other words, the birth of a cluster is a product of an uni- 
versal star formation process which operates on many scales 
of efficiency. 

In a recent work by GIO, it has been inferred that 
the present CFE (clusters formed in the last 10 Myr) is 
higher for a higher current SFR in the host, th e so-called F- 
log(EsFR) relation. ISilva- Villa fc LarsenI (|201l|) used a sam- 
ple of 5 nearby spiral galaxies to test the GIO relation using 
two different methods to estimate the CFRs. They observed 
that the recovered data points, in spite of the method used, 
scattered around the expected values and were impossible 
to reconcile with the Goddard et al. relation (see Figure [T]). 
Despite the discrepancy between the two results, we will 
include the r-log(EsFR) relation in the tests we will per- 
form for the BCGs. A discussion of the uncertainties affect- 
ing this relation will be presented in the next section. They 
may explain the disagreement between Goddard et al. and 
Silva- Villa fc Larsen results. 

To investigate whether BCGs follow these cluster-host 
relations we use the quantities listed in Table[T]and[2l Before 
we test the relation and compare our data to the ones pub- 
lished in the literature we discuss, briefly, the main source 
of uncertainties associated with the derived parameters and 
the used methods. 



3.1 Uncertainties 

3.1.1 The estimate o/ Esfr 

The estimate of Esfr depends on the SFR in the galaxy and 
the size of the host. The SFR is an averaged value of the star 
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Table 1. For each BCG, several quantities are listed in the table. From the left to the right column: name of the target; distance, absolute 
magnitude M^."^ '^'^ of the brightest young cluster in the galaxy; radius containing 80 % of the total B luminosity of the galaxy, Rgo%; 
mean SFR; surface density of SFR; CFR; and F. The SFR and CFR are estimated assuming a Salpeter IMF (see Section [XTTTJ ■ 



target 



distance^ A^^^Most 



R, 



■80% 



SFR 



CFR 



Mpc 



mag 



(Kpc) (Mgyr-l) (Moyr-lKpc-2) (Mgyr-i) 



(%) 



ESO 338-IG04 


37.5 


-15. 5=^ 


0.8 


3.2*^ 


1.55 


1.6 


50±10 


Haro 11 


82.3 


-16. 16=^ 


1.8 


22.0^= 


2.16 


11.2 


50 +\l^ 


ESO 185-IG13 


76.3 


-14.55^ 


1.98 


6.4<i 


0.52 


l.T^ 


26 ±5^^ 


Mrk930 


71.4 


-15.17<= 


1.7 


5.34° 


0.59 


1.33" 


25 ±10<= 


SBS 0335-052E* 


54 


-14.28 


0.66 


1.3 


0.95 


0.64 (0.14)s 


49 (10)s±15 



Ostlin et al 



Ostlinetal 



12003 ): 
1200] ): 



Adamo et al. 



Adamo et al 



AdamoetaL 
Reines et al.l 



2010a 1 
2011a ) 



12011^) : 
||2008| ') 



^ The va lues indicated between brackets are estimated for the 5 SSCs with ages ^ 10 Myr (see Figure[2land lAdamo et al.l 

I2OIOIJ) . However, in the analysis we use the values obtained including the whole cluster population. 
^ data from NED,|http://ned.ipac.caltech.edu/| 



Table 2. From the left to the right column: name of the target; 
extended radius, Rq, values in arcsec are in brackets; the specific 
luminosity, Tl, of the cluster population with respect the host in 
FUV (central wavelengths of the filter ~ 0.14 ^m), U, and B. 
In the case of SBS 0335-052E, we include between brackets the 
recovered Tl if only the 6 SSCs are considered. 



target 


Rg(Kpc) 


Tl(UV) 


Tl(U) 


Tl(B) 


ESO 338 


1.4 (7.9") 


- 


34.6 


20.9 


Haro 11 


4.9(12.5") 


20.3 


11.8 


16.4 


ESO 185 


4.9 (13.5") 


14.3 


9.8 


7.1 


Mrk930 


4.5 (12.5") 


19.2 


22.3 


12.1 


SBS 0335 


2.6 (10.0") 


29.4(21.3) 


32.8(24.0) 


38.4(26.3) 



formation happening in the system. Many tracers are used 
in the literature to estimate the SFR and they do not al- 
ways produce the same results. Moreover, in a galaxy there 
are regions which are quiescent and others starbursting. To 
estimate a meaningful Sspr value, it is hence necessary to 
determine the size of the region which is currently produc- 
ing stars. Since there is not a standard way to estimate the 
areas, the Esfr can be diluted or overestimated. When SFR 
or EspR are used to compare properties of sample of galaxies 
it is necessary to keep in mind that some of the scatter is 
due to the different methods applied to estimate these quan- 
tities. In the present analysis, we will compare our sample of 
BCGs with other literature data. For the BCGs, the SFRs 
are obtained from m easurements of the Ha fluxes (applying 
the iKennicuttlllQQa . law), and the sizes are estimated in a 
homogeneous way. 

Finally, as already pointed out by GIO, the Kennicutt 
law for SFR has been calibrated assuming a Salpeter IMF. 
However, when F is estimated, we compare CFR and SFR 
in the same host. Since in our star cluster analysis we have 
assumed a Kroupa IMF (2001), we have, for sake of consis- 
tency, applied a conversion factor in order to obtain cluster 
masses for a Salpeter IMF. The CFR and F values listed in 
Table [1] are for cluster masses derived assuming a Salpeter 
IMF. 



3.1.2 The area of the starburst 

The outskirts of BCGs are quite extended, while the star- 
burst regions are confined to the central areas of the galaxies. 
In the current analysis we use two different sizes (radii) : one 
which includes the starburst and outskirts {Kg, radius of the 
galaxy); the other to delimit the area where the star- forming 
regions are contained (R8o% . radius of the star-forming re- 
gion) . The Rg is the most extended radius we can infer from 
our data and is used to estimate the luminosity of the tar- 
gets, Lhost, in a few bands (see below). The R8o% is used for 
constraining the density of SFR in the galaxy, Esfr. The 
latter is a measure of the rate at which the star-forming re- 
gions are producing stars. Therefore, it is important to use 
a homogeneous method to estimate the area of these active 
regions when different galaxies are compared. 

To obtain an estimate of R8o% radius, we have looked 
at the fraction of galactic flux contained in growing aperture 
radii as function of the total flux contained inside Rg. A di- 
rect check on the frames shows that the R8o% encloses gen- 
erously the starburst regions in all the BCGs here studied. 
As a further check we have compared (Figure |3} the fraction 
of flux i n the B band and in Ha for 3 of the 5 BCGs (we 
refer to lOstlin et al.l |2009| for the reduction and analysis) . 
Ha emi ssion is usuall y considered a standard tracer of the 
SFR tKcnnicuttlll998h . On the other hand, the B band flux 
is dominated by the hght produced by young stars and can 
also be considered a reliable indicator of the location of the 
starburst regions. We compare the fraction of growing flux 
in the two tracers to check whether Rso% can be considered 
a good estimate of the size of the starburst regions in the 
galajcy (see Figure |3}. 

In the central regions, where the starburst dominate, 
the Ha luminosity distribution varies in the 3 galaxies. In 
ESO 338 (top panel, Figure|3}, we observe that the Ha lumi- 
nosity is less centrally concentrated than the B band lumi- 
nosity. However, at R8o%, we observe an opposite behaviour 
and roughly 90 % of the Ha flux is enclosed. This effect is 
caused by the presence of a very massive and young clus- 
ter in the centre of ESO 338 (|6stlin et al.ll2007^ . which has 
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Figure 3. The fraction of galactic luminosity contained inside 
increasing radii as function of the total luminosity at the most 
extended galactic radius, Rq. In the inset we show which target 
is analysed. The fraction of flux in 3 band, i(B) is showed as 
a black solid line, whereas the F(Ha) is the green (grey) dotted 
line. 



cleaned the surroun ding regions of the nebular gas (see im- 
age of the galaxy in lOstlin et al.ll2009l ). causing a dearth of 
Hq emission. In Haro 11 (center panel. Figure [3]), the situ- 
ation is inverted. The Hq luminosity is contained in a much 
smaller region than the B band one and, at R8o%, we de- 
tected almost the 100 % of Ha flux produce by the galaxy. 
Finally, in the case of SBS 0335-052E (bottom panel), we 
see that Ha and B band luminosities grow in a similar way. 
In all the cases within the uncertainties, Rgo%, esti- 
mated from the B band, incorporates or is in a fairly good 
agreement with Ha luminosity distribution. Therefore we 
use Rgo% to estimate the area of the star-forming regions in 
the galaxies and, thus, the Esfr. 



3.1.3 The distinction between hound clusters and 
associations 

The distinction between clusters and associations has not 
been applied to our BCGs sample, neither to the data used 
to obtain the three relations discussed in this paper. Two 
of these relations, the Tl(U)-Esfr and the r-log(EsFR), 
can be drastically affected. If the real number of clusters is 
overestimated (e.g. stars in unbound systems are counted as 



an evidence of an universal cluster formation efficiency, F of 
8 % if clusters form with a Schechter CMF. This relation 
also implies that galaxies with higher SFR form more mas- 
sive clusters due to a statistical sampling effect. However, 
if the cluster formation efficiency is constant for increasing 
SFR, the F and the Ti(U) should also be constant. There- 
fore, the positive relations of T^ (U) and F versus SFR could 
be caused by a contamination of unbound systems. There is 
no other evidence, which proves a universally constant clus- 
ter formation efficiency, nor is there an accessible way to 
estimate whether the T_l(U)-Esfr and the F-log(EsFR) re- 
lations still hold after a re-analysis including only bound 
objects. In this paper, we assume that clusters and asso- 
ciations are a product of the same star formation process, 
happening under different physical conditions and at differ- 
ent scales. This ass umption is also supported by the distri- 
bution presented in lGieles fc Portegies Zward (|20ld ). They 
see no clear distinction between associations and clusters 
during the first 10 Myr of their formation but a continu- 
ous distribution. Only when these systems age is there a 
clear distincti on with only clusters remaining tightly bound 
(see Fig. 2 bv lGieles fc Portegies Zwari]|2010l '). We refer to a 
cluster or association as a clustered str ucture, knowing that 
only bound clusters will survive longer (jBastian fc GoodwinI 
,2006i) and, eventually, become globular clusters. Since we are 
not able - with the current data - to disentangle these uncer- 
tainties, we will limit our analysis to compare the properties 
of the BCGs with these known relations and discuss possible 
implications, assuming that the estimated F and Tl(U) are, 
indeed, upper limits to the real values. 

3.1.4 Distance of the targets 

In general, the distance of the host system can affect the 
quality of the cluster analysis. For increasing distance, our 
ability to resolve single clusters diminishes rapidly. Blend- 
ing and crowding (clusters usually form in complexes) can 
significantly produce overestimates of the quantities we are 
interested in, Tl(U), m^"'^''''^"' , CFR and, thus, F. The clus- 
ter analysis of 3 of 5 BCGs are likely affected by blending, 
i.e., the most distant ones (distances are listed in Table [T]): 
Haro 11, ESO 185, Mrk930. The other two BCGs are close 
enough to resolve most of the clusters. Since derived clus- 
ter properties can be overestimated because of blending, we 
tried to look for any positive correlation between distance 
and F. In Figure 31 the distance of the hosts of the GIO 
sample and of the BCGs are plotted as function of their 
F. There is a trend between a high value of F and farther 
galaxies. However, number statistics are very poor and an 
observational bias affects this outcome. Low SFR systems 
are easier to observe locally, while their detection efficiency 
decreases rapidly at larger distances. On the other hand, 
high SFR hosts are easily observed at larger distances be- 
cause of their high luminosities. Another element that needs 
to be considered is that environments with SFRs compa- 
rable to the BCG ones are rare in the local volume (< 20 
Mpc). In this case, the position of the M83 starburst nu- 
cleus at about 5 Mpc and F > 20 % gives evidence that 
high star formation efficiency are also observed in systems 
where blending and crowding are not an issue. If we restrict 
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the comparison to the BCG sample, we do not see any dif- 
ference between the distant and nearby systems, supporting 
the idea that blending is not severe in either of the 3 most 
distant galaxies. 



3.2 The fraction of light in star clusters 

LROO introduced the specific luminosity for young star clus- 
ters defined as Tl ~ lOOLciustcrs/ihost- Tl gives an estimate 
of the fraction of the total galaxy light that is produced by 
stars in clustered regions. Using a sample of galaxies which 
includes quiescent spirals and star-forming dwarf systems, 
LROO found that the specific luminosity in the U band, 
Tl(U) and the SFR of the hosts were positively correlated. 
We present the same sample in Figure [5l including the BCG 
data points. We estimate Ldustors using only clusters for 
which a SED fit has been performed. The values of Tl(U), 
and for two other filters, FUV and B, are listed in Tabled 
The Tl(U) and the Esfr are higher in BCGs. In general, 
the trend suggests that the fraction of star formation hap- 
pening in clusters is important and increases as function of 
the SFR. 

To understand why we observe a much higher fraction 
of Tl(U) in BCGs than in spiral galaxies, we need to look at 
the different star formation histories in these two classes of 
hosts and not only at the SFRs. In general, we observe that 
the young star clusters in BCGs are preferentially clustered 
in clumps, in agreement with observations of high correla- 
tions in position as function of the age among young systems. 
Studies of spatial correlation among field stars, associations, 
and clusters show a higher clustering for youn ger samples, 
and a clear smoothing of the o lder structures (|Gieles et al.l 
I2OO8I : iBastian et al.ll2009l . l201ll ). The star formation in spi- 
ral systems has proceeded more or less constantly for a long 
lapse of time. Therefore cluster disruption has worked in 
favour of populating the stellar fields which form the bulk 
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Figure 5. Tj^(U)-Espr relation by LROO. The sample used by 
LROO is shown with black triangles. The BCGs arc added as pur- 
ple dots. 



of the optical luminosity in these systems. In BCGs, the 
starburst has been acting for rather short timescales (~ 40 
Myr) , suggesting that cluster disruption (mostly infant mor- 
tality) has not been effective. 

Looking at the values listed in Table O we see that 
the fraction of light produced by the star clusters increases 
at shorter wavelengths and contributes significantly to the 
UV and U luminosities of the BCGs. Therefore, it suggests 
that clustered regions in galaxies at redshift ~2-3, with 
metallicities similar to the BCGs, contribute to a consid- 
erable fraction of the UV-rest frame light. In studies of Ly- 
man brea k galaxy an alogs at redshift between ~0.1 and 0.2, 
[Ovcrzicr et al.l (|2008l ) observed super starburst compact re- 
gions which dominate the UV light of these targets. Simi- 
larl y compact clumps are also observed at redshift > 1 ga lax- 
ies JElmegreen et al.ll2009l : iForster Schreiber et al.ll201ll ). In 
our analysis, we find some evidence that those bright areas 
are probably unresolved star cluster knots recently formed 
and not dispersed yet by the interaction with the galactic 
environment, e .g. similarly to what is observed in BCGs an d 
in spiral arms (jBastian et al.ll2Q05l : lElmegreen et al.ll2006l ). 

3.3 Do BCGs follow the M^"«5''*'="*-SFR relation? 

The My"^ ^'^ -SFR relation was first noticed by L02, who 
suggested that the visual magnitude of the most luminous 
cluster and the total gala xy SFR where correlat ed. By means 
of numerical modelling, IWeidner et al.l (|2004 ) investigated 
the cluster physical conditions required to reproduce the re- 
lation, assuming that the brightest cluster would be also 
the most massive. They reproduced the observed relation 
under the assumption that a cluster population is formed 
within a time scale of ^ 10 Myr, following a CMF with in- 
dex steeper (~ —2.3) than is normally observed in cluste r 
populations (e.g. -2.0. IZhang fc Fall|[l999l : iBik et al.ll2003l '). 
BOS, releasing the condition of the brightest cluster also be- 
ing the most massive, observed that the relation showed less 
scatter if the brightest youngest (< 10 Myr) cluster was used 
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observed cluster luminosity function in the galaxy, assuming a 
disruption model (mass dependent disruption - MDD) to account 
for missing clusters. The BCGs are represented by dots (see inset). 



instead. Monte Carlo simulations of cluster populations with 
a CMF of power-law —2.0 or a Schechter CMF with power- 
law —2.0 and a characteristic mass of a few times IO^Mq 
where used to test different scenarios of cluster formation 
efficiency. Using the plot shown here in Figure [6l BOS ruled 
out a scenario where 100 % of the stars form in clusters with 
a CMF power-law of index —2.0. The L02 relation appeared 
to be in better agreement if clusters form with a Schechter 
CMF and only ~ 8 — 10 % of the stars reside in bound clus- 
ters (cluster formation efficiency, F ~ 0.1). We include the 
BCG sample in the My^'^ '^'^ -SFR, using the corresponding 
youngest brightest clusters. The 5 targets are located above 
the relation. The most nearby BCG in our sample, ESO 
338, is the one with the largest scatter. In agreement with 
the M^^f^^^^^'-SFR relation, a higher SFR in BCGs enables 
the formation of more massive (luminous) clusters than qui- 
escent spiral galaxies. Moreover, the position of the targets 
suggests that the F in BCGs is higher than the ~ 8 — 10 % 
found by B08. 

A higher CF E in BCGs has already been observed using 
the GIO relation jAdamo et al.ll2010al : lAdamo et al]|2011al : 
lAdamo et al.ll2011bn . In Figure [T] we show the relation in- 
cluding also SBS 0335 and ESO 338. The position of Haro 11 
and ESO 185 has slightly changed co mpared to the previou s 
publications (jAdamo et al.ll2010al and lAdamo et al.ll2011al ). 
because of the different IMF (Salpeter instead of Kroupa) 
and the different method used to estimate the area of nor- 
malization of the SFR (Section l3.H) . The data by Silva- Villa 
& Larsen are clearly scattered around the relation. BCGs 
follow the trend predicted by the relation, even if we notice 
that 3 of the 5 targets sit above the expected values. The 



mean value of the cluster formation efficiency in BCGs is 
FbcGs = 40 ± 10 %. 



4 CLUSTER FORMATION IN BCGS: A CLOSE 
LOOK AT HIGH CLUSTER FORMATION 
EFFICIENCIES 

We have used the cluster-host relations to investigate in 
more detail difference/analogies between BCGs and other 
star forming galaxies. Due the uncertainties affecting the 
data we consider the derived quantities as upper limits to 
the real values. 

One differe n ce of the sample of dwarf galaxies studied 
bv lBillett et al.l (|2002l ). luminous BCGs appear to be more 
efficient in forming star clusters. There are a few possibilities 
why this may be true. The Billet at al galaxies do not show 
any clear signature of merging events. Most likely cluster for- 
mation in these dwarf systems has been triggered by inter- 
nal instabilities, low shear, and inflows of gas. Moreover, the 
gravitational potential in these systems is not deep enough 
to retain the gas which is clearly observed in outflows. Lu- 
minous BCGs have recently experienced a merger event or 
accreted a considerable amo unt of gas (see lOstlin et al.ll200ll : 
lAdamo et ai]|2010al . r2011al lbh. This may have favoured the 
formation of hundreds of massive clusters. Moreover, these 
systems are likely one or two orders of magnitude more mas- 
sive than the ones in the Billett et al. sample, allowing them 
to retain the gas more efficiently. 

The My^'^ '^'^ -SFR relation, presented in the previous 
section, suggests that BCGs are, in the local universe, among 
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the systems which form very massive (bright) clusters. Al- 
though they follow the trend it is also clear that they are 
slightly offset from the prediction made by BOS of a con- 
stant cluster formation efficiency. The systematic scatter 
observed in the position of the BCGs in the M^''^ ^^'^'-SFR 
relation suggests that a higher CFE is operating in these 
systems, which are not dynamically relaxed, e.g., in a merg- 
ing/interaction phase. The trends observed in the other two 
relations, r-log(EsFR) and T_l(U)-Esfr, even if uncertain, 
support this scenario. A close look at the Goddard et al. rela- 
tion reveals an interesting point. In the diagram in Figure[7l 
two different groups can be delineated by the distributions 
of the data points. Using the values listed in Table 4 of GIO, 
we estimate that the mean CFE of the group with lower val- 
ues (the Silva- Villa & Larsen sample has been excluded) of 
EsFR is r = 8.7±4.3 %, close to the value found by BOS. The 
more efficient sample includes 2 targets from GIO (the star- 
burst nucleus of the spiral MS3 and the LIRG NCG 3256) 
and the BCGs. Their mean CFE is F = 35. 6± 10 %, roughly 
a factor of 4 higher. The two subsamples suggest that clus- 
ter formation efficiency is not constant at all scales of SFR. 
We see that in star-forming systems where star formation 
has proceeded more or less constantly without any signifi- 
cant burst, the CFE has a mean value of ~ 8 %. Starburst 
systems, on the other hand, are very active in producing 
star clusters. Possibly, the difference in CFEs refiects a dif- 
ferenc e in the conditio ns of the interstellar medium in the 
hosts (|Elmegreenll2008h . 

The universal cluster formation efficiency discussed by 
BOS could be valid in the local universe assuming that 
the interstellar medium density is lower than in high red- 
shift galaxies. However, it is difficult to extend its valid- 
ity to galaxies with an extreme environment (merging sys- 
tems). Numerical simulations have shown that GCs may 
have formed in strongly shocked media and high pres- 
sure fields, which have enhanced the gas compression and 
favoured the formation of more tightly bound str uctures 
JElmegreen fc Efremovl Il997l . iBournaud et al.1 120081 '). Such 
conditions are usually reached in galaxy mergers, where the 
very m assive young star c lusters a re observed (Anten nae 
system. [Mengel et al.ll2005l: Arp 2 20. IWilson et al.ll2006l : the 
Bird galaxv. IVaisanen et al.ll200Sl ). An increase of the CFE 
as function of redsh i ft has also been found in simulations by 
iMuratov fc GnedinI (|2010|). They observed that the cluster 
formation efficiency (e.g. the mass in cluster versus the total 
mass of the host) at redshift ~ 3 was much higher (about 
20 %) than then in the local universe. For this reason it is 
possible that the formation and survival of young and mas- 
sive clusters formed in extreme environments such as BCGs 
or LIRGs could help to trace the formation of the old GCs, 
if they have formed under similar conditions. 

However, only better data in terms of number statistics 
and quality (high resolution is required in order to distin- 
guish clusters from associations) are needed to reach more 
stringent conclusions. 



We have looked into the efficiency of the formation of 
massive star clusters. It is known that some properties of 
the cluster population and the mean star formation in the 
host system are correlated, suggesting that the formation of 
a cluster is not a local event but intrinsically connected to 
the mean properties of the galaxy. In general, we find that 
BCGs follow fairly well these relations, even if their SFRs 
and cluster properties are more extreme. 

We discuss possible uncertainties affecting our results. 
The relation which is least affected is the M^"^'^ ^''^'^'-SFR one, 
which suggests that the cluster formation efficiency is higher 
in BCGs than in quiescent spiral and dwarf starburst galax- 
ies. The same evidence is also suggested by the other two 
relations, r-log(EsFR) and Tl(U)-Esfr, which appear to be 
consistent with the general picture despite the uncertainties 
of the data. 

In particular, we observe that the inclusion of the BCG 
sample separates the r-log(EsFR) plane into two regions. 
Local spiral galaxies and dwarf starbursts as the Magellanic 
Clouds occupy the area around a mean CFE of 8 ± 5 % 
in agreement with the prediction made by BOS. The BCGs, 
together with a LIRG and a nuclear starburst region (in- 
cluded in the GIO sample), are in a region with a mean 
CFE of 35.6 ± 10 %. This indicates that the merger event 
has enhanced the cluster formation in these systems. 

We observe that the fraction of light produced by star 
clusters in BCGs increases at shorter wavelengths and con- 
tribute significantly to the UV and U luminosities. This sug- 
gests that clustered regions contribute a substantial fraction 
to the rest-frame UV light of z~l-3 galaxies with metallici- 
ties similar to the BCGs. In studies of Lyma n break galaxy 
analogs at redshift between ~0.1 and 0.2, lOverzier et al.l 
(|2008l ) observed super starburst compact regions which dom- 
inate the UV light of these targets. In our analysis, we find 
some evidence that those bright clumps are probably unre- 
solved star cluster knots, similarly to the ones observed in 
BCGs. 
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5 CONCLUSIONS 

To understand the role of BCGs in the galaxy formation 
scenario, we have tried to constrain how the star and cluster 
formation is operating in these systems. 
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